arXiv:l506.00042v2 [hep-ph] 15 Jul 2016 


A coupled-channel formalism for three-body final state interaction 

Peng Guo1’2.3.|^ 

^Physics Department, Indiana University, Bloomington, IN 47405, USA 
^Center For Exploration of Energy and Matter, Indiana University, Bloomington, IN 47408, USA. 
^Department of Physics and Engineering, California State University, Bakersfield, CA 93311, USA. 

(Dated: July 19, 2016) 

From dispersion relation approach, a formalism that describes final state interaction among three 
particles in a coupled-channel system is presented. Different representations of coupled-channel 
three-body formalism for spinless particles in both initial and final states are derived. 
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Introduction .—Hadron spectroscopy is one of impor- 
tant methods for studying non-perturbative QCD and 
gaining insights of hadron structures and decay mecha¬ 
nism. With high statistic data collected from facilities, 
such as BESIII, Jefferson Lab and Panda, data analysis 
becomes even more challenging than ever before, espe¬ 
cially, multiparticle dynamics may play the central pole 
in inelastic region, without proper consideration of con¬ 
straints from physics principles and multiparticle dynam¬ 
ics, amplitudes extracted from data may be misleading. 
Therefore, to understand phenomena precisely, theoret¬ 
ical description of decay amplitudes need to take into 
account all possible dynamics, and follow some basic 
physics principles, such as unitarity and analyticity. In 
the past, handling processes with multiple-particle final 
states has been mainly based on the isobar model mm, 
i.e. assuming a multiparticle decay proceeded through 
a series of quasi-two-body sequential decays. For exam¬ 
ple, a decay process of one particle (0) into three final 
states (1,2,3) is usually described by a sum of all possi¬ 
ble decay chains: 0 —>■ (12)3 -I- 1(23) -I- (31)2 —)• 123. For 
each individual decay chain, the amplitude is a product 
of kinematic factors, a coupling constant and a two-body 
amplitude that only depends on two-particle subenergy. 
Interaction among multiple final state particles has been 
ignored completely in the isobar model. 

Three-body correction to isobar model has been de¬ 
veloped since 1960 HE], which is based on subenergy 
dispersion relation approach by considering the unitarity 
and analyticity properties of amplitudes. In those disper¬ 
sive approaches [SHSmiHIIl, a decay amplitude is written 
as the sum of all possible decay chains. For each indi¬ 
vidual decay chain, the amplitude now is the product of 
kinematic factors, a subenergy dependent complex scalar 
function. This scalar function satisfies a coupled disper¬ 
sion relation equations, and the solutions of these equa¬ 
tions describe the rescattering effects among three parti¬ 
cles. In this approach, interaction among three particles 
is generated from pair-wise two-body interactions by ex¬ 
changing a particle between pairs. The unitarity and an- 
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alyticity are guaranteed naturally. However, all the pre¬ 
vious developments have not considered the contribution 
from inelastic channels yet. In reality, the subenergy in 
the most of hadron production processes usually is far be¬ 
yond the elastic region. Once inelastic channels open up, 
the interference between different channels may be im¬ 
portant nails]- In recent years, the demand for study¬ 
ing and including three-body effect has been increased 
significantly, such as, for excited baryon study at Jeffer¬ 
son Lab. The complication for establishing higher ex¬ 
cited baryon states in those studies are not only because 
most of those baryon states are produced from multiple- 
particle final states but also from the strongly coupled 
multiple channels in inelastic region. Similar situation 
may exist in incoming exotic mesons studies at Hall D, 
Jefferson Lab and ongoing excited charmonium studies at 
BES HI. Therefore, to disentangle all the coupled-channel 
effects from the multiple-particle final state interaction, 
a coupled-channel formalism for multiple-particle states 
is essential, some efforts based on effective theory formal¬ 
ism have been made along this line pDII22j . The goal of 
this work is to generalize dispersive three-body rescatter¬ 
ing formalism to include the channels in inelastic region. 
In this letter, the decay process of a spinless-particle to 
three spinless-particle is presented to demonstrate the 
basics of coupled-channel three-body formalism without 
complication of spin structure of particles. 

Basic representation of coupled-channel three-body 
formalism .—The decay of a spinless particle to three 
spinless particles is described by, 


(I(a)2(/3)3(7),out|0,in) = f(27r)^(5^( ^ p, - P)T“^3', 

2 ^ 1 , 2,3 

( 1 ) 

where {a, (3,^) stand for the species of final state parti¬ 
cles (1,2,3) respectively, and the four momenta of i-th 
final state particle and the parent particle are denoted 
by Pi and P respectively. The decay amplitude T°‘P'y are 
usually expressed as the sum of partial wave series in each 
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FIG. 1: A diagrammatic representation of discontinuity re¬ 
lations in Eq. the partial wave projection of Eq. gives 

Eq.0 . 


two-body subenergy-channel [3H9], 


+ PL(z;3^)Ff^^(s23)+PL(%a)Fr“'(s3l) , (2) 


( 7 “) / 


where the isospin couplings have been suppressed for 
simplification purpose only, the invariants are defined 
by Sij = (jpi -\-PjY s-nd three invariants are constrained 
by relation: si 2 -I- S 23 + S 31 = + Trip + rrii^ {M 

and m’s label parent and final state particle masses re¬ 
spectively). The total spin of two-particle subsystem is 
labeled by L. The cosine of polar angle of particle-1 in the 
rest frame of (l(a)2(/3)) system, Za^ = cosOap, is given 
by, 


812 ( 323 -S3i) +(ml - -m^) 

" 2 Mp^(si 2 ) 2 ^q^p{si 2 ) ’ ^ ’ 

where the momentum factors q and p are defined by 
, , a/[si2 - (rria - mpY] [si2 - (rUa + 

=-V5s-’ 

_ ^/[sl 2 -(M - m^Y] [si2 -(M + m^y] 
Pz(si 2 ) ^ ^ 

Similarly, the other z’s are given by cyclically permu- 
tating sub- and super-indices of Eqs. (§ and The 
dynamics of decay process are described by scalar func¬ 
tions E’s, which only depend on subenergy of isobar pair 
(sij) by assumption. 

Considering the analytic properties of decay amplitude 
the discontinuity crossing unitarity cut in suben¬ 
ergy, e.g. S 12 , then reads. 


Disci2T“^^(pi,P2,P3) 

= 1 V /" 

2(^0-^ (2^)32p;o (2^)32p'o 

X {2'n:f5'^(pi + P 2 - p'l - P 2 ) 

X {PlP2]p'lP2)T°‘'^'^ (Pl ,^2:^3) 

+ Y,i2L+l)PL(Zc.^)<j[^^\si2), (5) 


where the summation of (a'/3') run over all allowed two- 
body intermediate states for (V{a')2'{f3' j) pair, last term 
(jejiotes the contribution from the rest of inelastic 
channels. In current work, our discussion will be only 
limited to the three-body subspace of inelastic channels 
by choosing = 0. The partial wave expansion of 

two-body scattering amplitude AI in a coupled-channel 
system is given by 

Mal3^a'l5'(PlP2;p'lP2) 

= (IbTT) + ^)Pl(cos0) [iL(si2)](„^)_(„/,g/) , (6) 

L 


where 0 is the angle between incoming and outgo¬ 
ing particles of two-particle system. The matrix 
denotes coupled-channel partial wave scattering am¬ 
plitudes and it is normalized to Imt£^ = —p, where 
the non-vanishing elements of diagonal matrix p are 
given by A diagrammatic repre¬ 

sentation of discontinuity relations in Eq.([^ is shown 
in Fig. Commonly, the kinematical singularities 

are pulled out from decay amplitudes by defin¬ 
ing E^“^^(si 2 ) = F^l'^\si 2 )I [v^9a/3 ( 312 )^ 7 ( 512 )]^ 
[m [19], where possess only dynamical unitarity 

cuts by assumption. The discontinuity relations for 
scalar functions f\^^\si 2 ) are then derived from Eqs.([^ 
and (|^, 


Disci2Fl“^^ 

= E 

(«'/?') 

1 


( 312 ) 

K(3i2)6>(si 2 - 3fl)p(si2)](^^)_(^,^,) 
[ 7 ^ 90 / 3 ( 312 )^ 7 ( 312 )] ^ 

^ dZa'P'PL(Za'P>)T°‘ ^ {si 2 , 323, S 31 ), 


( 7 ) 


where the non-vanishing elements of diagonal matrix 
Sr are 3^^^ = (iria + mpY- The self-consistent integral 
equation for F^^\si 2 ) is constructed by dispersion re¬ 
lation. 


E, 


{aP) 


, Disci2E^“^^(s') 


(S12) = - 

TT 




ds‘ 


s' - Si 2 


( 8 ) 


where we have assumed that Disci 2 E^“^^(oo) = 0, so that 
no subtractions is needed. The angular projection in 
Eq.Q has to be analytically continued when disconti¬ 
nuity relation of Frs is plugged into Eq.([^, especially 
in the situation when the dispersion integral runs out of 
physical decay region, z’s are no longer defined on real 
axis between —1 and 1. The procedure of analytic con¬ 
tinuation has been given in mg. Similarly, sets of equa¬ 
tions for f\^"*\s 23 ) and Q^n be constructed 

in exactly the same approach, and together with Eq.(|^, 
they form a set of close coupled equations. The solu¬ 
tions of coupled-equation for F’s describe the three-body 
rescattering contribution from both elastic and inelastic 
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three-body channels. Eqs.Q and ([^ yield a basic rep¬ 
resentation of coupled-channel formalism for three-body 
final state interaction. The rescattering effect is produced 
by exchanging particle between isobar pairs, and the in¬ 
put of three-body equations are the two-body scattering 
amplitudes, which may be obtained from experimental 
measurements. 

Other representations of coupled-channel three-hody 
formalism .—Instead of solving Eqs. 0 and ([^, we may 
also consider other representations of three-body equa¬ 
tions, which demonstrate a explicit separation between 
rescattering contribution inside a pair and rescattering 
between pairs. 

As suggested in El, first of all, we may parametrize 
amplitudes by 


F^L^\s12) = [iL(Sl2)ffL(Sl2)](„^) , (9) 


where [f^] = [^L](a/ 3 ).(a'/ 3 ') / [9c/3ga'/3']^. In 

general, has both left-hand and right-hand singular¬ 
ities, i.e. Discti = 9{s — SL)ImtL q^^t*^9{s — Sfi)ptL, 
where sl labels branch points of left-hand singularities. 
Therefore, besides the unitarity cut, the vector has 
also left-hand singularities in order to keep F’s free off 
left-hand singularities, and discontinuity relations for 
thus read, 


Disci2 [5 l(si2)](„^) 

= - ^* l '^{ si 2)9 {sl - Sl2)ImfL(si2)gL(si2)](^^) 
+ ■:- 777^(512 — )Pap\S\2) 


[ 7 ^^ 7 ( 512 )] 

2L'F\ 


E 


/-I 


PL'izp^)F^^^'>\s23) + PL'MF^j^‘^’is3i) . (10) 


dZapPlj i^Zap) 


ha) , 


As illustrated in single channel case in E, when the dis¬ 
continuity relation for g is inserted into dispersion rela¬ 
tion, the integral equations in two variables are obtained: 
one variable is related to the angular projection; the other 
is associated to the dispersion integration. Fortunately, 
the Pasquier inversion technique [U HU [23] enable one to 
interchange the order of dispersive and angular integra¬ 
tions, and eventually write a single integral equations for 


5 L s. 


[5L(si2)](a/3) 
J°/3) 


1 

TT 


ds'.r 


s'l2 ~ ®12 


\t*L ^ImtL 5 L(si 2 )] 


(a/3) 


E 


p{M—ma)^ 


X \tL'{s 23 ) 9 L'{s 23 )] 


01 


dS3l S31) 

-00 

X [tL'(s3l)5L'(s3l)]. 


( 11 ) 


where the kernel functions /C„’s are defined in Eqs. (15 


161 , and /Cg’s do not depend on any dynamics but only 


on kinematic factors. Therefore, the ’universality’ prop¬ 
erties of ICg's allow one to compute them analytically, 
which is a great advantage for numerical evaluation of 
Eq.(ll). Similar equations for 31 ,( 523 ) and 3 l(s 3 i) are 
obtained by cyclic permutation of both sub- and super¬ 
indices in Eq.( 11). 

Next, we consider another representation of three-body 
equations by parameterization of 


p(a0) 

^ L 


(S 12 ) 


D{si2)Gl{Si2) 


( 12 ) 


where = t^Nff^ is denominator matrix functions of 
scattering amplitudes and has only right-hand singulari¬ 
ties by definition, and the left-hand singularities of tr are 
given by Ni matrix. and Nl are simply a coupled- 
channel generalization of standard N/D method PH dS] . 
In the single channel case, function niay be referred 
to as the Muskhelishvili-Omnes (MO) function PSI I27j . 
Thus, the vector G possess only right-hand singularities 
and the discontinuity relations for Gl read. 


Disci 2 [Gi,(si 2 )](„^) 

_ 9 q'/ 3 '(s 12 ) 

(a'/3') [^/^Pl{Sl2)] 


NI{si2)9{si2 - Sr)p{si2) 


(a/3).(a'/3') 


E 


2L' + 1 


j dZc^/p' Pij(^Zf^f p/') 


Pl'{zP'^)F^/'^\s 23) + PL'MF^r'\s3,) 


(13) 


Again, with the help of Pasquier inversion technique [S] 
ITU E5] . a single integral equations for Gl matrix are 
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obtained, 


= EE 

(a'P') U 


r(M-r. 


J (/3'7) 


dS 2 Z 

Di^} {s2z)GL'{S23,) 

/ (M —TTlo/)^ , , 

ds,i {s,2,s,,) 

-OO 

- (7«') 


D{s3i)GL’[ s^l) 


(14) 


where the kernel functions /Cg’s, together with kernel 
function /Cg’s defined in Eq.(ll), are given by 


rst, (s^i) foo 

e{s23) / (G') - 9{-S23) / (GO 

ds^,^{s23) •! st, (S2: 


■'g,G;LL' 

1 

TT 


9a'/3'(si2) [v^9/3'7('S23)Pa'(S23)] 


L' 


dSi2 

S'i2 — Sl2 


if 


Mp^{s[2) 




12^1 


(^2) 


(^' 12 ) 


2L' + 1 


J(a/3).(a'/3') 2 ^^' 7 ) ’ 

(15) 

^EEE^"“'^"'(^12,^3i) 

/-OO 1 J„/ 

0(S31)/ (GO -0(-a3i) / (GO 

,/(s3l) 


''g,G-,LL' 

1 

TT 


9a'/3'(s'l2) [v^97a'(s3l)4'/3'(s3l)] 


-W(33l) 

L' 


^\2 ~ -^12 


Mp^{s[2) 


12-^7 


(^ 12 ) 


if 


(S.G) 


(^ 2 ) 


2L' + 1 


J (a/3),(a'/3') 


PL{Za'P')PL'{zi^a')^ 


(16) 


where the matrix are given by -?f|f^(si 2 ) = ^ ^’^d 

d^L 0 ^ 12 ) = -^£(* 12 ) corresponding to gL and Gl respec¬ 
tively. The contour C is defined in Fig. 11 in [Mj, and 
the integration limits (the boundary of Dalitz plot), e.g. 
s^t,(s 23), are given by 


S^^(S23) = 


(ml 


-I- -I- -I- - S 23 

“ 2 ' 

- - ml) ^ 2Mpa{s23)qp-yis23) 

2S23 

( 17 ) 


Similar expression for s^c((®3i) are obtained by cyclically 
permutating indices in Eq. As we see in Eqs.([T^ 
and (16), the kernel functions Ac’s for7?L equations not 
only depend on dynamical functions Nl^s, but also has 
off-diagonal contributions from rescattering between elas¬ 
tic and inelastic channels due to non-diagonal matrix ■ 
As for gi^ equations, although, the kernel functions ICg’s 
are totally diagonal, the off-diagonal contributions ap¬ 
pear in the integral term over left hand cut (first term on 
the right-hand side of Eq.(ll)). Unlike ’universal’ kernel 
functions ICg, because of Nl dependence in kernel func¬ 
tions JCc, Ac’s now can only be computed by numerical 
intergration in complex plane. 


Finally, the integral equations for E, g and G pro¬ 
vide three equivalent representations of coupled-channel 
three-body formalism. As discussed in single channel 
three-body case in |14j . three different representations in 
principle yield the same result if II matrix is well-defined 
in complex plane. In practice, the information of tz, are 
usually only available in physical region on real axis, thus, 
different approximate methods for solving dispersion in¬ 
tegral equations are used. Therefore, the difference in 
solutions from different representation are expected de¬ 
pending on the approximations. In the single channel 
case |I4j . different approximate methods by restricting 
the integration ranges seem only change the overall nor¬ 
malization of solutions in physical region and barely alter 
the resonance properties, so the approximate solutions 
may be still justified. However, whether the conclusion 
still holds in coupled-channel case remains an open ques¬ 
tion. Nevertheless, single-integral-equation representa¬ 
tions for g and G are clearly easier to solve numerically 
and more suitable for event by event based data analysis. 


Summary .—In summary, we derived sets of integral 
equations for coupled-channel three-body final state in¬ 
teractions based on the dispersion approach, the for¬ 
malism is presented in three different representations in 
Eqs.([^, ([^ and ( [l4| . 
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